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Abstract 

Hybridization has played an important role in the evolutionary history of Canis 
species in eastern North America. Genetic evidence of coyote-dog hybridization 
based on mitochondrial DNA (mtDNA) is lacking compared to that based on 
autosomal markers. This discordance suggests dog introgression into coyotes 
has potentially been male biased, but this hypothesis has not been formally 
tested. Therefore, we investigated biparentally, maternally, and paternally inher- 
ited genetic markers in a sample of coyotes and dogs from southeastern 
Ontario to assess potential asymmetric dog introgression into coyotes. Analysis 
of autosomal microsatellite genotypes revealed minimal historical and contem- 
porary admixture between coyotes and dogs. We observed only mutually exclu- 
sive mtDNA haplotypes in coyotes and dogs, but we observed Y-chromosome 
haplotypes (Y-haplotypes) in both historical and contemporary coyotes that 
were also common in dogs. Species-specific Zfy intron sequences of Y-haplo- 
types shared between coyotes and dogs confirmed their homology and indicated 
a putative origin from dogs. We compared Y-haplotypes observed in coyotes, 
wolves, and dogs profiled in multiple studies, and observed that the Y-hapl- 
otypes shared between coyotes and dogs were either absent or rare in North 
American wolves, present in eastern coyotes, but absent in western coyotes. We 
suggest the eastern coyote has experienced asymmetric genetic introgression 
from dogs, resulting from predominantly historical hybridization with male 
dogs and subsequent backcrossing of hybrid offspring with coyotes. We discuss 
the temporal and spatial dynamics of coyote-dog hybridization and the condi- 
tions that may have facilitated the introgression of dog Y-chromosomes into 
coyotes. Our findings clarify the evolutionary history of the eastern coyote. 



Introduction 

The evolutionary histories and taxonomy of contempo- 
rary Canis species in North America are complex and 
controversial (see Nowak 2003; Chambers et al. 2012). 
There is agreement that the gray wolf (C. lupus) and 
coyote (C. latrans) are distinct species that evolved in 
Eurasia and North America, respectively (Nowak 1979; 
Kurten and Anderson 1980), and that the dog (C. famili- 
aris) was domesticated from the gray wolf (Wayne and 
Vila 2003). The eastern wolf (C. lycaon) is taxonomically 
controversial (e.g., vonHoldt et al. 2011 vs. Rutledge et al. 
2012a), but herein is considered an endemic North Amer- 
ican species based on consideration of genetic and nonge- 
netic data (Kyle et al. 2006; Mech 2011; Benson et al. 



2012; Chambers et al. 2012). Hybridization has clearly 
played an important role in the evolutionary history of 
Canis species in eastern North America (see Wayne and 
Vila 2003; Kyle et al. 2006; vonHoldt et al. 2011; Cham- 
bers et al. 2012). 

Hybridization is an evolutionary process that can occur 
naturally or because of anthropogenic influences. Regional 
habitat change that facilitates range expansion of one spe- 
cies into the range of another can result in hybridization, 
potentially producing novel gene combinations that can 
allow for rapid evolutionary change (Rhymer and 
Simberloff 1996). Such was the case for the coyote, which 
hybridized with the eastern wolf as it expanded into 
northeastern North America (Parker 1995; Wilson et al. 
2009; Kays et al. 2010). 
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Coyotes were historically confined to the western prai- 
ries and grasslands of North America (Young and Jackson 
1951), until the eradication of wolves and landscape 
changes associated with European colonization facilitated 
coyote expansion eastward during the last century (Parker 
1995). Coyotes differ noticeably in morphology and 
genetic composition between the western and eastern por- 
tion of the species' range. The larger size of eastern coyotes 
(Way 2007) has been attributed to either a phenotypic 
response to enhanced food supply (Thurber and Peterson 
1991) or a genotypic response to larger prey (Lariviere and 
Crete 1993), the latter of which was suggested to have been 
facilitated by genetic variation introduced through hybrid- 
ization with wolves (Kays et al. 2010). Eastern coyotes 
commonly exhibit introgressed wolf genes (Kays et al. 
2010; Way et al. 2010; Bozarth et al. 2011; Wilson 
et al. 2012), which are absent in western coyotes (Lehman 
et al. 1991; Hailer and Leonard 2008; Koblmiiller et al. 

2009) . The genetic evidence of wolf genes in eastern coy- 
otes corroborates earlier suggestions, based on cranial 
measurements of "New England Canis", that wolf genes 
had likely been introgressed into eastern coyotes (Law- 
rence and Bossert 1969; Mengel 1971). Furthermore, it 
was suggested that dog genes had likely been introgressed 
into eastern coyotes (Lawrence and Bossert 1969; but see 
Mengel 1971). Indeed, recent investigation of single- 
nucleotide polymorphisms (SNPs) has demonstrated dog 
admixture in eastern coyotes (vonHoldt et al. 2011), but 
it is intriguing that there is limited evidence of nonre- 
combining genes of dog origin in eastern coyotes, despite 
multiple studies investigating their genetic composition 
(Koblmiiller et al. 2009; Kays et al. 2010; Way et al. 2010; 
Wheeldon et al. 2010a; Bozarth et al. 2011). We are aware 
of only the observations of a single-dog mitochondrial 
DNA (mtDNA) haplotype in coyotes in the southeastern 
United States (Adams et al. 2003) and one putative dog- 
like mtDNA sequence in a coyote in Vermont (Kays et al. 

2010) . Thus, there is an apparent disparity in the genetic 
data pertaining to coyote-dog hybridization, which may be 
attributable to the focus on mtDNA in previous studies, or 
insufficient comparisons of genetic data across Canis spe- 
cies to infer the origins of introgression. It is plausible that 
genetic introgression from dogs to eastern coyotes may 
have occurred more extensively through paternal contribu- 
tions (e.g., Wheeldon and Patterson 2012), but this 
hypothesis has not been formally tested. Investigation of 
the Y-chromosome has proven valuable in assessing wolf- 
dog hybridization in Europe, revealing hybridization pat- 
terns that were not evident based on mtDNA (Iacolina 
et al. 2010; Godinho et al. 2011). Such could prove to be 
the case in assessing coyote-dog hybridization. 

Investigations of the canine Y-chromosome have typi- 
cally involved genotyping four to seven microsatellite loci 
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that are combined into Y-haplotypes (e.g., Sundqvist 
et al. 2001; Fain et al. 2010), but the recent evaluation of 
diagnostic SNPs in the last intron of the Zfy gene has 
facilitated the inference of the species-specific origins of 
these Y-haplotypes (Wilson et al. 2012). The combined 
investigation of Y-chromosome SNPs and microsatellite 
alleles, which are highly conserved (i.e., ~10~ 8 per site per 
generation in humans) and highly mutable (i.e., ~10~ 3 
per locus per generation in humans), respectively (see 
Hurles and Jobling 2001), facilitates the assessment of 
Y-haplotype diversity among species and the taxonomic 
origins of the genes. This approach is appropriate for the 
investigation of male-biased dog introgression into east- 
ern coyotes because of anticipated multispecies intro- 
gression (vonHoldt et al. 2011; Wilson et al. 2012). 

We generated autosomal microsatellite genotypes, 
mtDNA control region sequences, and Y-chromosome 
microsatellite genotypes and associated Zfy intron 
sequences for a sample of coyotes and dogs to investigate 
hybridization between these species in southeastern 
Ontario, Canada. The study region represents the putative 
geographic origin of the wolf-coyote hybridization that 
gave rise to the eastern coyote population that subse- 
quently colonized the maritime Canadian provinces and 
northeastern United States (Hilton 1978; Parker 1995; 
Wheeldon et al. 2010a), therefore the coyotes sampled 
herein should be representative of the greater eastern 
coyote population. Specifically, we employed mtDNA and 
Y-chromosome markers to test the hypothesis that eastern 
coyotes have experienced male-biased genetic introgres- 
sion from dogs, and we employed autosomal markers to 
assess levels of historical and contemporary gene flow 
between eastern coyotes and dogs. We note that our sam- 
pling of eastern coyotes is restricted to northeastern 
North America; reference to "eastern coyotes" with 
respect to our findings is intended to refer specifically 
northeastern coyotes. 

Material and Methods 

Samples and DNA extraction 

We obtained blood (n = 63), tissue (n = 273), and hair 
(« = 4) samples from 340 coyotes collected from south- 
eastern Ontario during 1974-1984 (n = 120) and 
2005-2010 (n = 220) (Fig. 1). The coyote samples were 
predominantly submitted by hunters and trappers or 
obtained from road kills. We generally restricted our sam- 
pling to south of 45° latitude, except for four samples 
collected east of 75° longitude and three collected west of 
81° longitude (Fig. 1), to reduce the potential of sampling 
eastern wolves in the vicinity of Algonquin Provincial 
Park. We obtained tissue (n = 16) and blood (n = 59) 
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Figure 1. Approximate locations of historical 
(1974-1984) and contemporary (2005-2010) 
coyote samples collected from southeastern 
Ontario, Canada; some samples may be 
overlapping. The black square in the inset map 
depicts the location of the study area in North 
America. Dog samples were collected in the 
Peterborough region. 




samples of 75 domestic dogs collected during 2006-2011, 
which were obtained predominantly from a local veteri- 
narian in the Peterborough region of Ontario (Fig. 1) 
and comprised mostly mixed-breed individuals. Breeds 
represented in the sampled dogs included Beagle, Border 
Collie, German Shepherd, Golden Retriever, Poodle, 
Husky, Labrador Retriever, Great Dane, Akita, Mastiff, 
Rottweiler, Malamute, Rhodesian Ridgeback, Spaniel, 
Terrier, and Hounds. 

We extracted DNA from samples using a DNeasy 
Blood and Tissue Kit (Qiagen, Mississauga, ON, Canada), 
and determined gender by amplification of Zfx/Sry primer 
pairs (P1-5EZ and P2-3EZ: Aasen and Medrano 1990; 
Y53-3C and Y53-3D: Fain and LeMay 1995) or Zfx/Zfy 
primers (LGL-331 and LGL-335: Shaw et al. 2003). 

Autosomal and Y-chromosome 
microsatellite genotyping 

For each sample we amplified 12 autosomal microsatellite 
loci in three multiplex reactions with published primers 
(cxx225, cxx2, cxxl23, cxx377, cxx250, cxx204, cxxl72, 
cxxl09, cxx253, cxx442, cxx410, and cxxl47: Ostrander 
et al. 1993, 1995) as in Wheeldon et al. (2010b). For male 
samples we amplified four Y-chromosome microsatellite 
loci with published primers (MS34A, MS34B, MS41A, 
and MS41B: Sundqvist et al. 2001), as in Wheeldon et al. 
(2010b). We purified amplified products through ethanol 
precipitation prior to genotyping on a MegaBACE 
1000 (GE Healthcare, Baie d'Urfe, QC, Canada) or an 
AB3730 (Applied Biosystems, Burlington, ON, Canada). 
We accounted for allele shifts between instruments with 



multiple control samples and scored alleles in Genemar- 
ker (vl.7, Softgenetics LLC, State College, PA). 

Mitochondrial DNA control region and Zfy 
intron amplification and sequencing 

We amplified a 343-347 base pair (bp) fragment of the 
control region of the mtDNA with published primers 
(AB13279: Pilgrim et al. 1998; AB13280: Wilson et al. 
2000) as in Wheeldon et al. (2010b). For some samples, 
we amplified the same fragment of the mtDNA control 
region with different published primers (ThrL, DL-Hcan: 
Leonard et al. 2002) under similar conditions to those in 
Wheeldon et al. (2010b). For a subsample of males we 
amplified a 658 bp fragment of the last intron of the Zfy 
gene with published primers (LGL-331: Shaw et al. 2003; 
Yint-2-335: Wilson et al. 2012) as in Wilson et al. (2012). 
We purified polymerase chain reaction (PCR) products 
using Exosap-IT (USB Corporation, Cleveland, OH), or 
Exonuclease I and Antarctic Phosphatase (New England 
BioLabs Inc., Ipswich, MA), prior to sequencing on a 
MegaBACE 1000 (GE Healthcare) or an AB3730 (Applied 
Biosystems). We edited and aligned sequences in Bioedit 
(v7.0.9, Hall 1999) or MEGA (v5, Tamura et al. 2011). 
Many of the coyote samples analyzed herein were previ- 
ously analyzed at the mtDNA control region by Wheel- 
don et al. (2010a). 

Genetic data analysis 

We obtained 411 autosomal microsatellite genotypes 
based on 12 loci for the coyote and dog samples; four 
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coyote samples (n = 3 historical; n = 1 contemporary) 
that amplified at less than six loci were excluded from 
subsequent analyses. Furthermore, we obtained genotypes 
of genetically assigned eastern wolves from Algonquin 
Provincial Park (n = 62: Rutledge et al. 2010) and gray 
wolves from northeastern Ontario (n = 62: Holloway 
2009; Rutledge et al. 2010; Wheeldon and Patterson 2012) 
for inclusion in autosomal data analyses to account for 
possible admixture from these groups in the sampled 
coyotes and to identify potential noncoyote migrants 
sampled in southeastern Ontario. We acknowledge the 
mixed ancestry of these wolf groups, but we refer to them 
as eastern wolves and gray wolves for simplicity. 

We analyzed the individual autosomal microsatellite 
genotypes in the Bayesian-clustering program Structure 
(v2.3, Pritchard et al. 2000; Hubisz et al. 2009) to assess 
admixture between coyotes and dogs. For all Structure 
analyses we inferred the parameter alpha and imple- 
mented the F-model (assumes correlated allele frequen- 
cies) and I-model (assumes independent allele 
frequencies) separately to compare results. We analyzed 
the historical and contemporary coyotes separately to 
avoid clustering problems associated with disparate sam- 
ple sizes among groups. Based on prior findings (Rutledge 
et al. 2010; Benson et al. 2012) and accounting for the 
inclusion of dogs, we anticipated that K = 4 would be 
optimal given the data. To confirm this for the data set 
that included contemporary coyotes, eastern wolves, gray 
wolves, and dogs, we ran the admixture model of Struc- 
ture for K = 1-7 with five repetitions of 10 iterations fol- 
lowing a burn-in period of 10 5 iterations for each K. We 
calculated the mean posterior probability (In P[D]) for 
each K by averaging across the five runs and confirmed 
that K = 4 was optimal for the data set based on quanti- 
tative criteria (Fig. SI; In P[D], Pritchard et al. 2000; AK, 
Evanno et al. 2005) and consideration of the biological 
significance of clusters. Additionally, we performed a fac- 
torial correspondence analysis on the entire data set of 
individual autosomal microsatellite genotypes (n = 535) 
as implemented in GENETIX (Belkhir et al. 2004) and 
observed clustering patterns generally concordant with the 
results from Structure (Fig. S2). Notably, there was com- 
plete overlap of the historical and contemporary coyote 
samples (Fig. S2), indicating that ^-determination in 
Structure was justifiably unnecessary for the data set that 
included historical coyotes. Subsequently, for both data 
sets, we ran the admixture model of Structure 10 times at 
K = 4 for 10 6 iterations following a burn-in period of 10 5 
iterations and collected information on the 90% probabil- 
ity intervals of individual assignments (ANCEST- 
DIST =1). We obtained individual admixture 
proportions (Q-values) from the run with the highest 
posterior probability and lowest variance; we considered 
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individuals to be admixed, if Q < 0.8 (e.g., Vaha and 
Primmer 2006). We compared individual assignments 
between the F-model and I-model and observed general 
concordance (Fig. S3). Specifically, we observed five cases 
of an individual being assigned as admixed under the 
F-model but not admixed under the I-model model, and 
five cases of the reverse scenario (Fig. S3). We suggest 
neither model was optimal for our data set considering 
the variable evolutionary relationships between Canis spe- 
cies, therefore we averaged Q-values between the F-model 
and I-model; hereafter references to individual assign- 
ments are based on the averaged Q-values, unless stated 
otherwise. 

We generated 223-228 bp mtDNA sequences and dis- 
tinguished those evolved in gray wolves and dogs (i.e., 
Old World origin) from those evolved in coyotes and 
eastern wolves (i.e., New World origin) based on a diag- 
nostic indel (Pilgrim et al. 1998; Wilson et al. 2000). 

We generated 400 bp Zfy intron sequences and 
observed three previously described sequences based on 
two variable sites: Yint-1 of coyote origin; Yint-2 of gray 
wolf (or dog) origin; and Yint-4 of eastern wolf origin 
(Wilson et al. 2012). We did not observe Yint-3 of coyote 
origin (Wilson et al. 2012) in our sample, which probably 
reflects a founder effect whereby eastward colonizing 
coyotes carried the Yint-1 sequence. 

We generated Y-haplotypes based on the Y-chromo- 
some microsatellite genotypes. We constructed a median- 
joining network (Bandelt et al. 1999) based on combined 
Y-chromosome microsatellite and Zfy intron data using 
the program Network (v4.6.1.0; available at www.fluxus 
-engineering.com/sharenet.htm) with default settings 
(e = 0). The Y-chromosome microsatellite loci were 
weighted equally and the Zfy intron sequence variation 
was weighted twice as high as microsatellite loci. 

We compared the Y-haplotypes that we observed in 
coyotes and dogs in this study with those of previously 
analyzed wolves, coyotes, and dogs to investigate Y-haplo- 
type sharing among these species for the purpose of infer- 
ring introgression. We defined western and eastern 
coyotes as occurring west and east of the Mississippi 
River, respectively, and Great Lakes coyotes as occurring 
in the western Great Lakes states and western Ontario 
(Wheeldon et al. 2010b). We defined Great Lakes wolves 
as occurring in Manitoba, northern Ontario, southern 
Quebec, and the western Great Lakes states (Wheeldon 
2009). Allele sizes for this study were calibrated with 
those in Rutledge et al. (2010), Wheeldon et al. (2010b), 
Wheeldon and Patterson (2012), and Wilson et al. 
(2012), the last of which calibrated allele sizes with those 
in Hailer and Leonard (2008). Furthermore, allele sizes 
were standardized among Musiani et al. (2007), Hailer 
and Leonard (2008), Sundqvist et al. (2001, 2006) (allele 
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sizes obtained from the authors), and Fain et al. (2010). 1 
Thus, allele sizes were standardized for comparisons 
across studies. We generated unique Y-haplotype identifi- 
ers (Yhn) for each Y-chromosome microsatellite genotype 
because existing identifiers were not the same for match- 
ing Y-haplo types across studies (Table SI). We performed 
a reduced-median analysis (r = 2; Bandelt et al. 1995) of 
the Y-haplotypes (loci weighted equally; frequency >1) 
and then constructed a median-joining network (e = 10) 
with maximum parsimony postprocessing (Polzin and 
Daneshmand 2003). 

Results 

Autosomal microsatellite data 

We observed in coyotes varying levels of autosomal 
admixture from dogs, eastern wolves, and/or gray wolves 
(Fig. 2). One contemporary coyote was identified as an 
eastern wolf migrant 2 (Q EW > 0.8) and another was iden- 
tified as eastern wolf x dog admixed with negligible 
coyote assignment (Fig. 2); both individuals were omitted 
from further consideration. Approximately 94% of histor- 
ical (n = 110) and 89% of contemporary (n = 194) coy- 
otes had individual assignments with Qcoy ^ 0.8, and 
~3% of historical (n = 3) and contemporary (« = 6) coy- 
otes had individual assignments with Qdog > 0.2 (Fig. 2). 
Approximately 2% of historical (n = 2) and 4% of 
contemporary (« = 9) coyotes had individual assignments 
with Q EW > 0.2, and no historical (n = 0) and 1% of 
contemporary (n = 3) coyotes had individual assignments 
with Qgw > 0.2 (Fig. 2). All the historical and contempo- 
rary coyotes with Qcoy > 0.8 had 90% probability inter- 
vals for Qcoy that did not overlap zero (both models), 
whereas only two historical coyotes (both models) and 
one contemporary coyote (F-model only) with 
Qdog > 0.2 had 90% probability intervals for Q D og that 
did not overlap zero (Table S2). Autosomal genotypes of 
several coyotes that had individual assignments with 
Qdog > 0.2 exhibited 1-2 alleles that were relatively 
common in dogs but rare in coyotes, corroborating that 
admixture had occurred. 

Mitochondrial DNA sequence data 

We observed 11 and 15 mutually exclusive mtDNA hapl- 
otypes in coyotes and dogs, respectively (Table S3). The 

'The published allele sizes and frequencies of multiple Y-haplo- 
types were erroneous; we obtained the corrected allele sizes and 
frequencies from the authors. 

2 Male individual omitted for Y-haplotype frequency determina- 
tion in coyotes. 



mtDNA haplotypes observed in dogs were exclusively of 
Old World origin (i.e., dog or gray wolf) and those 
observed in coyotes were exclusively of New World origin 
(i.e., coyote or eastern wolf). 

Y-chromosome data 

We obtained Y-chromosome microsatellite genotypes for 
218 male samples and observed 20 Y-haplotypes. Specifi- 
cally, we observed 15 and 10 Y-haplotypes in coyotes and 
dogs, respectively; five were shared between coyotes and 
dogs (Table 1). We obtained Zfy intron sequences for all 
Y-haplotypes and found that those observed exclusively in 
coyotes were associated with Yint-1 (n = 6), Yint-2 
(n = 1), or Yint-4 (n = 3), and those observed exclusively 
in dogs were associated exclusively with Yint-2 (Table 1). 
The Y-haplotypes shared between coyotes and dogs were 
associated exclusively with Yint-2 (Table 1). The Y-haplo- 
type network revealed a clear separation between the 
coyote-specific and dog-specific Y-haplotypes (Fig. 3); we 
note that Y-haplotype HT has been observed in dogs 
(Table 2), therefore we consider it shared regardless of it 
not being observed in dogs profiled in this study. The 
Y-haplotypes shared between coyotes and dogs clustered 
with the dog-specific Y-haplotypes (Fig. 3), therefore we 
consider them as belonging to the dog Y-haplogroup. 
Accordingly, there was a minimum of three microsatellite 
mutational steps and a Zfy intron substitution between 
the coyote and dog Y-haplogroups (Fig. 3). 

Comparison of the Y-haplotypes observed in coyotes 
and dogs in this study with those observed in wolves, coy- 
otes, and dogs in previous studies revealed Y-haplotype 
sharing variously between and among species (Table 2). 
We observed Y-haplotypes in coyotes that were shared with 
dogs (« = 7), eastern wolves (n = 3), Great Lakes wolves 
(n = 5), North American gray wolves (n = 3), and Eur- 
asian gray wolves (n = 2); some of these Y-haplotypes were 
shared jointly with two wolf groups, or with wolves and 
dogs. Eastern coyotes shared more Y-haplotypes with dogs 
(n = 7) than did Great Lakes coyotes (n = 3) or western 
coyotes (n = 0). We also observed Y-haplotypes that were 
shared exclusively between and among wolf groups, or 
between wolves and dogs. 

Discussion 
Autosomal admixture 

We observed a low proportion (3%) of coyotes exhibiting 
autosomal admixture from dogs (Fig. 2), and the propor- 
tion was the same during the historical and contemporary 
sampling periods, suggesting that coyote-dog hybridiza- 
tion has been minimal during the past 35-40 years in 
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■ Coyote Eastern Wolf ■ Gray Wolf ■ Dog 




1974-1984 2005-2010 



Figure 2. Structure assignments of historical (1974-1984) and contemporary (2005-2010) coyotes from southeastern Ontario genotyped at 12 
autosomal microsatellite loci. Each partitioned vertical bar represents an individual's proportional membership to the K = 4 populations analyzed 
in Structure; Q-values represent the average of the F-model and l-model results. 

Table 1. Y-chromosome haplotypes and associated Zfy intron sequences observed in historical (1974-1984) and contemporary (2005-2010) 
coyotes from southeastern Ontario and dogs. Haplotype frequencies are provided and the number of samples that had the Zfy intron sequenced 
for each haplotype is indicated in parentheses. Haplotypes are grouped by Zfy intron and sorted based on the allele at locus MS41A. 



Haplotype 


MS34A 


MS34B 


MS41A 


MS41B 


Zfy intron 


Historical coyotes 


Contemporary coyotes 


Dogs 


CD 


172 


178 


214 


210 




28 (0) 


47 (6) 




cw 


172 


178 


214 


208 






1 (D 




AJ 


172 


180 


212 


214 






2(2) 




CR 


172 


178 


212 


216 




4(2) 


2(2) 




GO 


176 


180 


212 


220 




1 (D 


4(2) 




GP 


176 


180 


212 


222 




7(1) 


11 (2) 




BB 


170 


182 


212 


226 


4 


1 (D 


2 (1) 




BQ 


170 


182 


212 


218 


4 




1 (D 




AA 


172 


180 


212 


212 


4 


13 (1) 


19(1) 




FE 


174 


178 


208 


216 


2 


1 (D 




1 (D 


FF 


174 


178 


208 


222 


2 




7(3) 


7(3) 


FG 


174 


178 


208 


224 


2 




1 (D 


1 (D 


FL 


174 


178 


208 


218 


2 


2(2) 


6(2) 


8(3) 


FS 


174 


178 


208 


226 


2 






10 (3) 


FT 


174 


178 


208 


220 


2 






6(3) 


HC 


174 


176 


208 


214 


2 






1 (D 


HG 


174 


176 


208 


224 


2 






1 (D 


HT 


174 


176 


208 


220 


2 




13 (3) 




JT 


174 


172 


208 


220 


2 


1 (D 


5(3) 


3 (3) 


FZ 


174 


178 


196 


220 


2 






1 (D 



southeastern Ontario. We observed a low proportion 
(4%) of coyotes exhibiting autosomal admixture from 
eastern wolves or gray wolves (Fig. 2), concordant with 
our sampling of coyotes outside of primary wolf range in 
Ontario; some of the coyotes that exhibited autosomal 
admixture from wolves could have been dispersing indi- 
viduals from central Ontario where contemporary wolf- 
coyote hybridization occurs (Rutledge et al. 2010; Benson 
et al. 2012). The proportion of coyotes exhibiting autoso- 
mal admixture from eastern wolves or gray wolves was 
higher in the contemporary group (5%) than the histori- 
cal group (2%), which likely reflects geographical differ- 
ences in sampling between time periods, rather than a 
change in the frequency of wolf-coyote hybridization, 
because there were many contemporary but few historical 
samples collected from locations near where wolves 
occurred (Fig. 1; Benson et al. 2012). 



Introgression of nonrecombining markers 

We observed mtDNA haplotypes of coyote or eastern 
wolf origin, but not those of dog or gray wolf origin, in 
our sample of coyotes (Table S3), which is concordant 
with previous studies (e.g., Koblmuller et al. 2009; Way 
et al. 2010; Bozarth et al. 2011; but see Adams et al. 
2003; Kays et al. 2010). Notably, mtDNA haplotype C13 
was only observed in one individual, which assigned pre- 
dominantly as eastern wolf in Structure (Table S2), 
consistent with the lack of observation of this coyote-clus- 
tering sequence in coyotes and it being part of the eastern 
wolf lineage (Wheeldon and White 2009; Fain et al. 
2010). 

We observed Y-haplotype sharing between dogs and 
coyotes (Table 1, 2), which we attribute to introgression 
from dogs into coyotes. We reject homoplasy as a poten- 
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Figure 3. Median-joining network of Y-chromosome haplotypes observed in coyotes from southeastern Ontario and dogs. Haplotypes are 
composites of Y-chromosome microsatellite genotypes and Zfy intron sequences and are represented by nodes that have sizes proportional to 
their frequency: gray, coyote; black, dog. Median vectors are represented by white nodes. Red text on branch lengths indicates the number of 
mutational steps between nodes: N, nucleotide substitution; R, microsatellite repeat. 



tial cause of the observed Y-haplotype sharing between 
dogs and coyotes because the Zfy intron sequences of the 
shared Y-haplotypes confirm their homology and putative 
origin from dogs (Tables 1, 2; Wilson et al. 2012). Fur- 
thermore, multiple mutational steps would need to have 
occurred in coyote Y-chromosomes for homoplasy to 
have generated the Y-haplotypes shared with dogs 
(Fig. 3), which seems unlikely. We reject the hypothesis 
that the observed Y-haplotype sharing between coyotes 
and dogs is attributable to a shared common ancestor 
given that they are present in eastern coyotes but appar- 
ently absent in western coyotes (Table 2), the latter of 
which represent the source of eastward colonizers that 
evolved into the former; this implies that coyotes experi- 
enced genetic introgression during eastward range expan- 
sion. Similarly, the eastern-specific presence in coyotes of 
the Y-halpotypes shared between coyotes and gray wolves 
or eastern wolves implies that they occur in coyotes due 
to introgression (Wilson et al. 2012); probably derived via 
hybridization with eastern wolves (Rutledge et al. 2010). 

The Y-haplotype sharing that we observed between 
wolves and dogs (Table 2) could variously be the result of 
shared ancestry, homoplasy (e.g., Fig. S4; Sundqvist et al. 
2006; Hailer and Leonard 2008), historical hybridization 
(Vila et al. 1997), or contemporary hybridization (e.g., 
Godinho et al. 2011; Hindrikson et al. 2012), but formal 
assessment of these alternatives is beyond the scope of 
this study. Despite sharing multiple Y-haplotypes, wolves 



and dogs also exhibited multiple mutually exclusive 
Y-haplotypes (Table 2). Notably, of the seven Y-haplo- 
types shared between coyotes and dogs, none were 
observed in North American gray wolves and eastern 
wolves, one was observed in Great Lakes wolves, and two 
were observed in Eurasian gray wolves (Table 2), the last 
of which are not a plausible source of introgression for 
coyotes because of geographical considerations. The con- 
temporary haplotype diversities of North American 
wolves (i.e., gray wolves, Great Lakes wolves, and eastern 
wolves) generally do not reflect their historical haplotype 
diversities because of genetic drift associated with drastic 
changes in population sizes (Leonard et al. 2005; Leonard 
and Wayne 2008). Accordingly, more of the Y-haplotypes 
shared between coyotes and dogs may have occurred in 
North American wolves historically, such that potentially 
these Y-haplotypes could have been introgressed into coy- 
otes from wolves historically. However, gray wolves are 
aggressive toward and kill coyotes (e.g., Thurber et al. 
1992; Berger and Gese 2007; Merkle et al. 2009), and 
genetic data indicate that these species do not interbreed 
in the wild (e.g., Pilgrim et al. 1998; Koblmiiller et al. 
2009); we are not aware of any successful interbreeding 
between these species in captivity. Similarly, Great Lakes 
wolves are aggressive toward and kill coyotes (see Mech 
2011) and rarely hybridize with them (e.g., Fain et al. 
2010; Wheeldon et al. 2011b). Eastern wolves readily 
hybridize with coyotes (Rutledge et al. 2010) and there- 
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fore potentially could have mediated gene flow from dogs 
to coyotes historically, but there is no evidence of dog 
introgression in eastern wolves sampled to date (Table 2; 
Rutledge et al. 2012b; Benson et al. 2012b). Coyotes and 
dogs have been interbred in captivity (Gier 1968; Silver 
and Silver 1969; Mengel 1971); direct mating between 
dogs (i.e., domesticated gray wolves) and coyotes in the 
wild may be facilitated because they are not ecological 
competitors and because of similarities in size and 
appearance of specific dog breeds with respect to coyotes. 
We suggest that the Y-haplotypes shared between coyotes 
and dogs were introgressed into coyotes from dogs, not 
wolves. Genetic profiling of historical samples may or 
may not reveal Y-haplotype sharing trends not observed 
in contemporary samples and analysis of additional mark- 
ers could resolve differences between apparently shared 
Y-haplotypes (e.g., Brown et al. 2011), providing either 
support for our hypothesis or alternative hypotheses. 

Temporal and spatial dynamics of 
introgression 

Considering the contrasting levels of admixture detected 
based on the autosomal and Y-chromosome markers ana- 
lyzed in this study, we suggest the present diversity of int- 
rogressed dog genes in the greater eastern coyote 
population originated predominantly from historical 
rather than contemporary hybridization between these 
species. This is supported by our observation that 19% 
(n = 31 of 166) of male coyotes with Qcoy 5 0.8 exhib- 
ited a dog Y-haplotype (Table S2), but only 3% of coy- 
otes exhibited autosomal admixture from dogs 
(Qdog > 0.2, Fig. 2). Specifically, 6% (n = 3 of 54) and 
25% (n = 28 of 112) of historical and contemporary male 
coyotes with Qcoy ^ 0.8, respectively, exhibited a dog 
Y-haplotype (Table S2). This indicates that the frequency 
of introgressed dog Y-chromosomes in coyotes in south- 
eastern Ontario increased between 1974-1984 and 2005- 
2010 despite no apparent change in autosomal admixture 
from dogs (Fig. 2). Notably, five of the six male coyotes 
with Qdog > 0.2 exhibited a dog Y-haplotype (Table S2: 
n = 1 historical; n = 4 contemporary); these individuals 
putatively represent Fl hybrids or backcrosses. The dog 
admixture detected in northeastern coyotes based on 
SNPs is estimated to have originated ~30 years ago 
(vonHoldt et al. 2011), but our Y-chromosome data dem- 
onstrated dog introgression in highly assigned coyotes 
sampled from southeastern Ontario >30 years ago (Table 
S2), indicating that hybridization had occurred earlier. 
This discrepancy may reflect geographical differences in 
sampling of eastern coyotes between studies and the like- 
lihood that coyotes continued to interbreed with dogs at 
the periphery of their range as they expanded eastward 



(Hilton 1978; Parker 1995). We suggest that the intro- 
gressed dog Y-haplotypes observed in eastern coyotes 
originated from direct coyote-dog hybridization that 
occurred primarily historically, and concomitantly with 
wolf-coyote hybridization in southeastern Ontario, as 
coyotes expanded eastward during the last century, 
although coyote-dog hybridization may also have 
occurred during eastward expansion across the western 
Great Lakes states. Colonizing coyotes expanding at the 
eastern periphery of the species' range would presumably 
have occurred at relatively low densities and may have 
been subject to an Allee effect (Allee 1931), encountering 
free-roaming dogs near human settlements in greater fre- 
quency than conspecifics, resulting in hybridization. This 
scenario has been documented for wolf-dog hybridiza- 
tion, whereby hybridization occurred in areas of low wolf 
density (Andersone et al. 2002; Munoz-Fuentes et al. 

2010) or was restricted to peripheral and recently expanded 
wolf populations (Godinho et al. 2011). The nonrecombin- 
ing dog Y-chromosomes were probably introgressed into 
the expanding coyote population at low levels and then 
amplified in frequency by logistic demographic growth that 
presumably occurred as coyote density increased and con- 
specific breeding predominated (Currat et al. 2008); this 
could explain the common occurrence of dog Y-haplotypes 
in contemporary eastern coyotes despite minimal evidence 
of autosomal admixture from dogs. The eastern coyote 
may exhibit varying levels of contemporary autosomal 
admixture from dogs across its range. 

Introgressed dog genes are common in coyotes in 
northeastern North America and those occurring east of 
the Mississippi River and south of the Great Lakes 
(Table 2; vonHoldt et al. 2011; Wilson et al. 2012), which 
may reflect the aforementioned susceptibility of an 
expanding population to hybridization. Morphological 
evidence of coyote-dog hybridization in western North 
America has been reported (Bee and Hall 1951; Mahan 
et al. 1978; Freeman and Shaw 1979), but available 
genetic data do not support the introgression of dog 
genes into coyotes in that region (Table 2; vonHoldt et al. 

2011) , suggesting that backcrossing of such hybrids with 
coyotes has been rare there. The dog (and wolf) Y-haplo- 
types observed in Great Lakes coyotes could be the result 
of in situ hybridization or gene flow from westward dis- 
persing eastern coyotes. Additional Y-chromosome profil- 
ing could further elucidate the geographic extent of 
coyotes that exhibit introgressed dog Y-haplotypes. 

Conditions facilitating introgression 

Coyotes and dogs interbred in captivity can produce fer- 
tile offspring (Gier 1968; Silver and Silver 1969; Mengel 
1971), and wild coyote-dog hybrids have been reported 
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(Gipson et al. 1974; Mahan et al. 1978; Freeman and 
Shaw 1979). However, it has been postulated that dog 
genes were unlikely to become introgressed into the 
coyote gene pool because a phase shift in the breeding 
season of coyote-dog hybrids makes backcrossing with 
coyotes unlikely (Gier 1968; Silver and Silver 1969; Men- 
gel 1971). Additionally, it has been suggested that lack of 
parental care from male dogs or coyote-dog hybrids 
would reduce hybrid offspring survival (Silver and Silver 
1969; Mengel 1971). Our findings indicate that despite 
these potential inhibiting factors, some hybrid male off- 
spring from female coyote x male dog crosses must have 
successfully backcrossed with coyotes, as noted for off- 
spring from a female dog x male coyote crossing (Adams 
et al. 2003). The backcrossing of wolf-dog hybrids with 
wolves has also been reported (Godinho et al. 2011). The 
introgression of dog Y-chromosomes into coyotes may 
have been possible because of variability in the breeding 
season of male coyote-dog hybrids that could have facili- 
tated backcrossing with female coyotes (Gipson et al. 
1975). Furthermore, pack associates may have helped pro- 
vision for hybrid offspring (see Andelt 1985) in the 
absence of male parental care, which does not preclude 
mateless female coyotes potentially raising pups to inde- 
pendence (Sacks and Neale 2001). 

Causes of asymmetric introgression 

The predominantly asymmetric nature of coyote-dog 
hybridization parallels that of wolf-dog hybridization, 
whereby hybridization typically involves a female of the 
wild species and a male of the domestic species (see God- 
inho et al. 2011). Evidence of female dog x male wolf or 
coyote crossings is relatively rare (Adams et al. 2003; 
Muhoz-Fuentes et al. 2010; Hindrikson et al. 2012), 
which could reflect a lack of recruitment of such hybrid 
offspring into wild populations. Female dogs, excepting 
feral ones, are unlikely to give birth in the wild, therefore 
their hybrid offspring would be unlikely to backcross with 
the wild species and presumably would not be sampled in 
wildlife genetic surveys (Hindrikson et al. 2012). How- 
ever, the different breeding cycles of male and female 
dogs may facilitate asymmetric hybridization between 
wild canids and dogs, considering that male dogs can 
breed year round and therefore are capable of mating 
with female wild canids in heat, but female dogs only 
come into heat twice a year in variable seasons and there- 
fore only a fraction of them are capable of mating with 
reproductively active male wild canids (Mengel 1971). 
Asymmetric introgression from dogs into coyotes may 
result partly from the potentially reduced fecundity of 
female coyote-dog hybrids (Mengel 1971; Gipson et al. 
1975) and the potentially decreased fitness of coyote-dog 



hybrids with dog mtDNA relative to those with coyote 
(or eastern wolf) mtDNA. Notably, Haldane's Rule (Hal- 
dane 1922) does not apply for coyote-dog hybrids. 

Diagnostic value of locus MS41A 

The value of locus MS41A as a diagnostic marker for distin- 
guishing Y-haplotypes of gray wolf (and dog) origin from 
those of coyote (and eastern wolf) origin has been alluded to 
previously (Hailer and Leonard 2008), and recently Zfy 
intron sequences have confirmed species-specific origins of 
Y-haplotypes (Wilson et al. 2012). Based on our findings 
(Table 1, 2) and those of Wilson et al. (2012), Y-haplotypes 
with allele 208 at locus MS41A are associated exclusively 3 
with Yint-2 of gray wolf or dog origin, and those with alleles 
212, 214, or 216 at locus MS41A are associated with Yint-1 
(alleles 212, 214, and 216) or Yint-3 of coyote origin (alleles 
212 and 214), or Yint-4 of eastern wolf origin (allele 212). 
Additionally, a rare Y-haplotype with allele 196 at locus 
MS41A, which was observed exclusively in dogs, was associ- 
ated with Yint-2 (Table 1). TheZfyintron(s) associated with 
Y-haplotypes with allele 210 or 218 at locus MS41A remain 
to be determined. Analysis of additional Y-chromosome 
SNPs is required for potentially differentiating Y-haplotypes 
that are apparently shared between gray wolves and dogs 
(e.g., Brown et al. 201 1) because the Zfy intron fragment we 
sequenced lacks resolution in this respect. 

Conclusions 

Our investigation of the canine Y-chromosome has further 
clarified the evolutionary history of the eastern coyote. Spe- 
cifically, Y-haplotypes observed in eastern coyotes that were 
previously attributed to gray wolf introgression have now 
been predominantly attributed to dog introgression, 
although some gray wolf introgression is still evident 
(Table 2). We observed three incomplete Y-chromosome 
genotypes in coyotes in this study that appeared to be of gray 
wolf origin (Table S2), therefore we potentially underesti- 
mated the level of gray wolf introgression in eastern coyotes. 
Although allele sizes from Koblmuller et al. (2009) were 
unavailable for comparison, we question whether the eight 
"wolf clade" Y-haplotypes observed exclusively in eastern 
coyotes in that study were of wolf origin; we suggest the pos- 
sibility that those Y-haplotypes were of dog origin. 

The evidence presented herein of introgressed dog and 
wolf genes in eastern coyotes is corroborated by investiga- 
tions of skull morphology (Lawrence and Bossert 1969; 
McGinnis 1979; Kays et al. 2010) and coat color (Ander- 



3 A Y-haplotype with allele 208 at locus MS41A observed in a 
Texas coyote was associated with Yint-3 and was suggested to 
result from a rare homoplasy (Wilson et al. 2012). 



3016 



© 201 3 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd. 



T. J. Wheeldon ef al. 



Dog Y-Chromosome Introgression into Coyotes 



son et al. 2009; Brockerville et al. 2013). The introgres- 
sion of wolf genes has contributed to the large body size 
and skull dimensions of eastern coyotes, putatively 
making them more effective predators of deer (Kays et al. 
2010), but the effect, if any, of introgressed dog genes on 
the evolution of eastern coyotes remains unknown. We 
speculate that behavioral traits of dogs would likely be 
selected against in coyotes because loss of wildness would 
presumably lead to increased mortality risk from humans. 
We are not aware of any evidence to suggest that the 
introgression of wolf or dog genes has resulted in more 
aggressive behavior of eastern coyotes toward humans. 

This study and others (e.g., Iacolina et al. 2010; 
Wheeldon et al. 2012) highlight the importance of using 
multiple genetic markers when assessing hybridization 
between species. The standard approach for studies inves- 
tigating hybridization should involve the assessment of 
maternally, paternally, and biparentally inherited genetic 
markers, and if possible should incorporate morphologi- 
cal data (e.g., Benson et al. 2012; Wheeldon and Patter- 
son 2012). 

Acknowledgments 

We thank B. Noble for providing dog samples, the many 
hunters and trappers from southeastern Ontario for 
contributing coyote samples, and NSERC for funding. 

Data Accessibility 

DNA sequences: Genbank accessions FJ687618-FJ687620 
(Zfy intron) and Table S3 (mtDNA). Microsatellite data 
deposited in the Dryad Repository: http://dx.doi.org/ 
10.5061/dryad.72s5k. 

Author Contributions 

TJW conducted laboratory work and data analyses and 
wrote the manuscript. All coauthors contributed to inter- 
pretation of results and writing of manuscript. 

Conflict of Interest 

None declared. 
References 

Aasen, E., and J. F. Medrano. 1990. Amplification of the Zfy 
and Zfx genes for sex identification in humans, cattle, sheep 
and goats. Biotechnology 8:1279-1281. 

Adams, J. R., J. A. Leonard, and L. P. Waits. 2003. 

Widespread occurrence of a domestic dog mitochondrial 
DNA haplotype in southeastern US coyotes. Mol. Ecol. 
12:541-546. 



2013 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd. 



Allee, W. C. 1931. Animal aggregations: a study in general 
sociology. University of Chicago Press, Chicago, IL. 

Andelt, W. F. 1985. Behavioral ecology of coyotes in south 
Texas. Wildl. Monogr. 94:3-45. 

Anderson, T. M., B. M. vonHoldt, S. I. Candille, M. Musiani, 

C. Greco, D. R. Stahler et al. 2009. Molecular and 
evolutionary history of melanism in North America gray 
wolves. Science 323:1339-1343. 

Andersone, Z., V. Lucchini, E. Randi, and J. Ozolins. 2002. 
Hybridisation between wolves and dogs in Latvia as 
documented using mitochondrial and microsatellite DNA 
markers. Mamm. Biol. 67:79-90. 

Bandelt, H.-J., P. Forster, B. C. Sykes, and M. B. Richards. 
1995. Mitochondrial portraits of human populations using 
median networks. Genetics 141:743-753. 

Bandelt, H.-J., P. Forster, and A. Rohl. 1999. Median-joining 
networks for inferring intraspecific phylogenies. Mol. Biol. 
Evol. 16:37^8. 

Bee, J. W., and E. R. Hall. 1951. An instance of coyote-dog 
hybridization. Trans. Kansas Acad. Sci. 54:73-76. 

Belkhir, K, P. Borsa, L. Chikhi, N. Raufaste, and 

F. Bonhomme 2004 GENETIX 4.05, logiciel sous windows 
TM pour la genetique des populations, laboratoire genome, 
populations, interactions, CNRS UMR 5171, Universite de 
Montpellier II, France. 

Benson, J. F., B. R. Patterson, and T. J. Wheeldon. 2012. 
Spatial genetic and morphologic structure of wolves and 
coyotes in relation to environmental heterogeneity in a 
Canis hybrid zone. Mol. Ecol. 21:5934-5954. 

Berger, K. M., and E. M. Gese. 2007. Does interference 
competition with wolves limit the distribution and 
abundance of coyotes? J. Anim. Ecol. 76:1075-1085. 

Bozarth, C. A., F. Hailer, L. L. Rockwood, C. W. Edwards, and 
J. E. Maldonado. 2011. Coyote colonization of northern 
Virginia and admixture with Great Lakes wolves. 
J. Mammal. 92:1070-1080. 

Brockerville, R. M., M. J. McGrath, B. L. Pilgrim, and 
H. D. Marshall. 2013. Sequence analysis of three 
pigmentation genes in the Newfoundland population of 
Canis latrans links the Golden Retriever Mclr variant to 
white coat color in coyotes. Mamm. Genome. doi:10.1007/ 
s00335-012-9443-x. 

Brown, S. K, N. C. Pederson, S. Jafarishorijeh, 

D. L. Bannasch, K. D. Ahrens, J.-T. Wu, et al. 2011. 
Phylogenetic distinctiveness of Middle Eastern and 
Southeast Asian village dog Y chromosomes illuminates 
dog origins. PLoS One 6:e28496. 

Chambers, S. M., S. R. Fain, B. Fazio, and M. Amaral. 2012. 
An account of the taxonomy of North American wolves 
from morphological and genetic analyses. North Am. Fauna 
77:1-67. 

Currat, M., M. Ruedi, R. J. Petit, and L. Excoffier. 2008. The 
hidden side of invasions: massive introgression by local 
genes. Evolution 62:1908-1920. 

3017 



Dog Y-Chromosome Introgression into Coyotes 

Evanno, G., S. Regnaut, and J. Goudet. 2005. Detecting the 

number of clusters of individuals using the software 

STRUCTURE: a simulation study. Mol. Ecol. 14:2611-2620. 
Fain, S. R., and J. P. LeMay. 1995. Gender identification of 

humans and mammalian wildlife species from PCR amplified 

sex linked genes. Proc. Am. Acad. Forensic Sci. 1:34. 
Fain, S. R., D. J. Straughan, and B. F. Taylor. 2010. Genetic 

outcomes of wolf recovery in the western Great Lakes states. 

Conserv. Genet. 11:1747-1765. 
Freeman, R. C. and J. H. Shaw. 1979. Hybridization in Cants 

(Canidae) in Oklahoma. Southwest. Nat. 24:485^99. 
Gier, H. T. 1968. Coyotes in Kansas. Bull, agric. exp. stat., 

Kansas State University 393:1-118. 
Gipson, P. S., J. A. Sealander, and J. E. Dunn. 1974. The 

taxonomic status of wild Canis in Arkansas. Syst. Zool. 

23:1-11. 

Gipson, P. S., I. K. Gipson, and J. A. Sealander. 1975. 

Reproductive biology of wild Canis (Canidae) in Arkansas. 

J. Mammal. 56:605-612. 
Godinho, R., L. Llaneza, J. C. Blanco, S. Lopes, F. Alvares, 

E. J. Garcia, et al. 2011. Genetic evidence for multiple events 

of hybridization between wolves and domestic dogs in the 

Iberian Peninsula. Mol. Ecol. 20:5154-5166. 
Hailer, F., and J. A. Leonard. 2008. Hybridization among three 

native North American Canis species in a region of natural 

sympatry. PLoS One 3:e3333. 
Haldane, J. B. S. 1922. Sex ratio and unisexual sterility in 

hybrid animals. J. Genet. 12:101-109. 
Hall, T. A. 1999. Bioedit: a user-friendly biological sequence 

alignment editor and analysis program for Windows 95/98/ 

NT. Nucleic Acids Symp. Ser. 41:95-98. 
Hilton, H. 1978. Systematics and ecology of the eastern coyote. 

Pp. 209-228 in M. Bekoff, ed. Coyotes: biology, behavior, 

and management. The Blackburn Press, Caldwell, NJ. 
Hindrikson, M., P. Mannil, J. Ozolins, A. Krzywinski, and 

U. Saarma. 2012. Bucking the trend in wolf-dog 

hybridization: first evidence from Europe of hybridization 

between female dogs and male wolves. PLoS One 7:e46465. 
von Holdt, B. M., J. P. Pollinger, D. A. Earl, J. C. Knowles, 

A. R. Boyko, H. Parker et al. 2011. A genome-wide 

perspective on the evolutionary history of enigmatic 

wolf-like canids. Genome Res. 21:1294-1305. 
Holloway, J. 2009. Pp. 81 Size dependent resource use of a 

hybrid Wolf (C. lycaon x C. lupus) population in Northeast 

Ontario. [Masters thesis], Trent University, Peterborough. 
Hubisz, M. J., D. Falush, M. Stephens, and J. K. Pritchard. 

2009. Inferring weak population structure with the 

assistance of sample group information. Mol. Ecol. Resour. 

9:1322-1332. 

Hurles, M. E., and M. A. Jobling. 2001. Haploid chromosomes 
in molecular ecology: lessons from the human Y. Mol. Ecol. 
10:1599-1613. 

Iacolina, L., M. Scandura, A. Gazzola, N. Cappai, C. Capitani, 
L. Mattioli, et al. 2010. Y-chromosome microsatellite 



3018 



T. J. Wheeldon et al. 

variation in Italian wolves: a contribution to the study of 

wolf-dog hybridization patterns. Mamm. Biol. 75:341-347. 
Kays, R. K., A. Curtis, and J. J. Kirchman. 2010. Rapid 

adaptive evolution of Northeastern coyotes via hybridization 

with wolves. Biol. Lett. 6:89-93. 
Koblmiiller, S., M. Nord, R. K. Wayne, and J. A. Leonard. 

2009. Origin and status of the Great Lakes Wolf. Mol. Ecol. 

18:2313-2326. 

Kurten, B., and E. Anderson. 1980. Pleistocene mammals of 

North America. Columbia Univ. Press, NY. 
Kyle, C. J., A. R. Johnson, B. R. Patterson, P. J. Wilson, 

K. Shami, S. K. Grewal, et al. 2006. Genetic nature of 

eastern wolves: past, present and future. Conserv. Genet. 

7:273-287. 

Lariviere, S., and M. Crete. 1993. The size of eastern coyotes 
(Canis latrans): a comment. J. Mammal. 74:1072-1074. 

Lawrence, B., and W. H. Bossert. 1969. The cranial evidence 
for hybridization in New England Canis. Brevoria 330:1-13. 

Lehman, N., A. Eisenhawer, K. Hansen, L. D. Mech, 

R. O. Peterson, P. J. P. Gogan, et al. 1991. Introgression of 
coyote mitochondrial DNA into sympatric North American 
gray wolf populations. Evolution 45:104-119. 

Leonard, J. A., and R. K. Wayne. 2008. Native Great Lakes 
wolves were not restored. Biol. Lett. 4:95-98. 

Leonard, J. A., R. K. Wayne, J. Wheeler, R. Valadez, 

S. Guillen, and C. Vila. 2002. Ancient DNA evidence for old 
World origin of new World dogs. Science 298:1613-1616. 

Leonard, J. A., C. Vila, and R. K. Wayne. 2005. Legacy lost: 
genetic variability and population size of extirpated US grey 
wolves (Canis lupus). Mol. Ecol. 14:9-17. 

Mahan, B. R., P. S. Gipson, and R. M. Case. 1978. 

Characteristics and distribution of coyote x dog hybrids 
collected in Nebraska. Am. Midi. Nat. 100:408-415. 

McGinnis, H. J 1979. Pp 227 Pennsylvania coyotes and their 
relationship to other wild Canis populations in the Great 
Lakes region and the northeastern United States. [Masters 
thesis], Pennsylvania State University, State College. 

Mech, L. D. 2011. Non-genetic data supporting genetic 
evidence for the eastern wolf. Northeast. Nat. 18:521—526. 

Mengel, R. M. 1971. A study of dog-coyote hybrids and 

implications concerning hybridization in Canis. J. Mammal. 
52:316-336. 

Merkle, J. A., D. R. Stahler, and D. W. Smith. 2009. 

Interference competition between gray wolves and coyotes 
in Yellowstone National Park. Can. J. Zool. 87:56-63. 

Muhoz-Fuentes, V., C. T. Darimont, P. C. Paquet, and 
J. A. Leonard. 2010. The genetic legacy of extirpation and 
re-colonization in Vancouver Island wolves. Conserv. Genet. 
11:547-556. 

Musiani, M., J. A. Leonard, H. D. Cluff, C. C. Gates, 
S. Mariani, P. C. Paquet, et al. 2007. Differentiation of 
tundra/taiga and boreal coniferous forest wolves: genetics, 
coat colour and association with migratory caribou. Mol. 
Ecol. 16:4149-4170. 



© 201 3 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd. 



T. J. Wheeldon ef al. 



Dog Y-Chromosome Introgression into Coyotes 



Nowak, R. M. 1979. North American quaternary canis. 

Monogr. Mus. Nat. Hist. 6:1-154. 
Nowak, R. M. 2003. Wolf evolution and taxonomy. Pp. 

239-258 in L. D. Mech and L. Boitani, eds. Wolves: 

behavior, ecology, and conservation. University of Chicago 

Press, Chicago, IL. 
Ostrander, E. A., G. F. Sprague, and J. Rine. 1993. 

Identification and characterization of dinucleotide repeat 

(CA)n markers for genetic mapping in dog. Genomics 

16:207-213. 

Ostrander, E. A., F. A. Mapa, M. Yee, and J. Rine. 1995. One 
hundred and one new simple sequence repeat-based markers 
for the canine genome. Mamm. Genome 6:192-195. 

Parker, G. 1995. Eastern coyote: the story of its success. 
Nimbus Publishing, Halifax, NS. 

Pilgrim, K. L., D. K. Boyd, and S. H. Forbes. 1998. Testing for 
wolf-coyote hybridization in the Rocky Mountains using 
mitochondrial DNA. J. Wildl. Manag. 62:683-689. 

Polzin, T., and S. V. Daneshmand. 2003. On Steiner trees and 
minimum spanning trees in hypergraphs. Oper. Res. Lett. 
31:12-20. 

Pritchard, J. K., M. Stephens, and P. Donnelly. 2000. Inference 

of population structure from multilocus genotype data. 

Genetics 155:945-959. 
Rhymer, J. M., and D. Simberloff. 1996. Extinction by 

hybridization and introgression. Annu. Rev. Ecol. Syst. 

27:83-109. 

Rutledge, L. Y., C. J. Garroway, K. M. Loveless, and 
B. R. Patterson. 2010. Genetic differentiation of eastern 
wolves in Algonquin Park despite bridging gene flow 
between coyotes and grey wolves. Heredity 
105:520-531. 

Rutledge, L. Y., P. J. Wilson, C. F. C. Kliitsch, B. R. Patterson, 
and B. N. White. 2012a. Conservation genomics in 
perspective: a holistic approach to understanding Canis 
evolution in North America. Biol. Conserv. 155:186-192. 

Rutledge, L. Y., B. N. White, J. R. Row, and B. R. Patterson. 
2012b. Intense harvesting of eastern wolves facilitated 
hybridization with coyotes. Ecol. Evol. 2:19-33. 

Sacks, B. N, and J. C. C. Neale. 2001. Does paternal care of 
pups benefit breeding female coyotes? Southwest. Nat. 
46:121-126. 

Shaw, C, P. Wilson, and B. White. 2003. A reliable molecular 
method of gender determination for mammals. J. Mammal. 
84:123-128. 

Silver, H., and W. T. Silver. 1969. Growth and behavior of the 
coyote-like canid of northern New England with 
observations on canid hybrids. Wildl. Monogr. 17:3-41. 

Sundqvist, A.-K., H. Ellegren, H. Olivier, and C. Vila. 2001. 
Y chromosome haplotyping in Scandinavian wolves (Canis 
lupus) based on microsatellite markers. Mol. Ecol. 
10:1959-1966. 

Sundqvist, A.-K., S. Bjornerfeldt, J. A. Leonard, F. Hailer, A. 
Hedhammar, H. Ellegren, et al. 2006. Unequal contribution 



of sexes in the origin of dog breeds. Genetics 
172:1121-1128. 

Tamura, K., D. Peterson, N. Peterson, G. Stecher, M. Nei, and 
S. Kumar. 2011. MEGA5: molecular evolutionary genetics 
analysis using maximum likelihood, evolutionary distance, 
and maximum parsimony methods. Mol. Biol. Evol. 
28:2731-2739. 

Thurber, J. M., and R. O. Peterson. 1991. Changes in body 
size associated with range expansion in the coyote (Canis 
latrans). J. Mammal. 72:750-755. 

Thurber, J. M., R. O. Peterson, J. D. Woolington, and 
J. A. Vucetich. 1992. Coyote coexistence with wolves 
on the Kenai Peninsula, Alaska. Can. J. Zool. 70: 
2494-2498. 

Vaha, J. P., and C. R. Primmer. 2006. Efficiency of 

model-based Bayesian methods for detecting hybrid 

individuals under different hybridization scenarios and with 

different numbers of loci. Mol. Ecol. 15:63-72. 
Vila, C, P. Savolainen, J. E. Maldonado, I. R. Amorim, 

J. E. Rice, R. L. Honeycutt, et al. 1997. Multiple and ancient 

origins of the domestic dog. Science 276:1687-1689. 
Way, J. G. 2007. A comparison of body mass of Canis latrans 

(coyotes) between eastern and western North America. 

Northeast. Nat. 14:111-124. 
Way, J. G., L. Rutledge, T. Wheeldon, and B. N. White. 2010. 

Genetic characterization of eastern "coyotes" in eastern 

Massachusetts. Northeast. Nat. 17:189-204. 
Wayne, R. K., and C. Vila. 2003. Molecular genetics studies of 

Wolves. Pp. 218-238 in L. D. Mech and L. Boitani, eds. 

Wolves: behavior, ecology, and conservation. University of 

Chicago Press, Chicago, IL. 
Wheeldon, T. 2009. Pp.120 Genetic characterization of Canis 

populations in the Western Great Lakes Region. [Masters 

thesis], Trent University, Peterborough. 
Wheeldon, T. J., and B. R. Patterson. 2012. Genetic and 

morphological differentiation of wolves (Canis lupus) and 

coyotes (Canis latrans) in northeastern Ontario. Can. J. 

Zool. 90:1221-1230. 
Wheeldon, T., and B. White. 2009. Genetic analysis of 

historical western Great Lakes region wolf samples reveals 

early Canis lupus/ lycaon hybridization. Biol. Lett. 

5:101-104. 

Wheeldon, T., B. Patterson, and B. White. 2010a. Colonization 
history and ancestry of northeastern coyotes. Biol. Lett. 
6:246-247. 

Wheeldon, T., B. Patterson, and B. White. 2010b. Sympatric 
wolf and coyote populations of the western Great Lakes 
region are reproductively isolated. Mol. Ecol. 
19:4428-4440. 

Wheeldon, T., B. Patterson, and D. Beyer. 2012. Coyotes in 
wolves' clothing. Am. Midi. Nat. 167:416-420. 

Wilson, P. J., S. Grewal, I. D. Lawford, J. N. M. Heal, A. G. 
Granacki, D. Pennock et al. 2000. DNA profiles of 
the eastern Canadian wolf and the red wolf provide 



2013 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd. 



3019 



Dog Y-Chromosome Introgression into Coyotes 



T. J. Wheeldon ef al. 



evidence for a common evolutionary history independent 
of the gray wolf. Can. J. Zool. 78:2156-2166. 

Wilson, P. J., S. K. Grewal, F. F. Mallory, and B. N. White. 
2009. Genetic characterization of hybrid wolves across 
Ontario. J. Hered. 100:S80-S89. 

Wilson, P. J., L. Y. Rutledge, T. J. Wheeldon, B. R. 

Patterson, and B. N. White. 2012. Y-chromosome evidence 
supports widespread signatures of three-species Canis 
hybridization in eastern North America. Ecol. Evol. 
2:2325-2332. 

Young, S. P., and H. H. T. Jackson. 1951. The clever coyote. 
The Stackpole Company, Harrisburg, PA. 

Supporting Information 

Additional Supporting Information may be found in the 
online version of this article: 

Figure SI. Plots of quantitative JC-determination criteria, 
In P(D) and AK, for K = 1-7 Structure analyses of con- 
temporary (2005-2010) coyotes from southeastern 
Ontario, eastern wolves from Algonquin Provincial Park, 
gray wolves from northeastern Ontario, and dogs geno- 
typed at 12 autosomal microsatellite loci. The F-model 
and I-model were implemented separately. 
Figure S2. Factorial correspondence analysis of historical 
(1974-1984) and contemporary (2005-2010) coyotes from 
southeastern Ontario, eastern wolves from Algonquin 
Provincial Park, gray wolves from northeastern Ontario, 
and dogs genotyped at 12 autosomal microsatellite loci. 
Figure S3. Structure assignments of contemporary 
(a = F-model; b = I-model) and historical (c = F-model; 
d = I-model) coyotes from southeastern Ontario, eastern 
wolves from Algonquin Provincial Park, gray wolves from 



northeastern Ontario, and dogs genotyped at 12 autoso- 
mal microsatellite loci. Each partitioned vertical bar 
represents an individual's proportional membership to 
the K = 4 populations analyzed in Structure. 
Figure S4. Median-joining network of Y-chromosome 
haplotypes observed in coyotes, wolves, and dogs. Haplo- 
types are based on Y-chromosome microsatellite geno- 
types and are represented by variably colored nodes: WC, 
yellow; GLC, green; EC, orange; EW, red; GLW, purple; 
NA GW, blue; EA GW, gray; and DD, black (abbrevia- 
tions defined in Table 2). Node sizes are not proportional 
to haplotype frequency. Median vectors are represented 
by white nodes. Branch lengths are not proportional to 
the number of mutational steps between nodes; branch 
lengths represent one microsatellite repeat (R), except 
where indicated otherwise. 

Table SI. Y-chromosome haplotype identifiers provided 
in different studies: ID1, Rutledge et al. 2010; Wheeldon 
et al. 2010b; Wheeldon and Patterson 2012; and Wilson 
et al. 2012; ID2, Hailer and Leonard 2008; ID3, Fain et al. 
2010; ID4, Sundqvist et al. 2001; ID5, Sundqvist et al. 
2006. 

Table S2. Sample information, genetic data, and Structure 
assignments (K = 4) and associated 90% probability inter- 
vals (PI) of historical (1974-1984) and contemporary 
(2005-2010) coyotes from southeastern Ontario genotyped 
at 12 autosomal microsatellite loci. Structure Q-values: red, 
averaged; yellow, F-model; blue, I-model. Abbreviations: 
COY, coyote; EW, eastern wolf; GW, gray wolf. 
Table S3. Frequencies and Genbank accessions of mito- 
chondrial DNA haplotypes observed in historical 
(1974-1984) and contemporary (2005-2010) coyotes from 
southeastern Ontario and dogs. 
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